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Відновлення гравіпотенціалу за модулем його градієнта в задачі Алексідзе 

Пропоную вашій увазі результати дослідження нової нелінійної граничної задачі гравіметрії, яка змінює усталені по- 
няття про аномалії сили тяжіння і матиме за вдалої реалізації широке коло застосувань (слайо Г). 

Для вирішення задач вивчення фігури, внутрішньої будови, глибинної геодинаміки Землі слід знати розподіл значень 
потенціалу сили тяжіння, значень градієнта цього потенціалу або значень модуля цього градієнта. Цей розподіл можна от- 
римати з фундаментальних властивостей функцій, що описують гравітаційну взаємодію, 1 вимірів. Створена теорія віднов- 
лення потенціалу в глобальній області за слідами потенціалу або слідами його похідних на границі області. Але здійснен- 
ню глобальних побудов заважають наступні проблеми і труднощі (слайд 2): 

е відсутні точні граничні дані для реалізації класичної схеми відновлення в глобальній області потенціалу як розв'язку однієї з гра- 
ничних задач Діріхле, Неймана чи змішаної (Стокса-Молоденського) для рівняння Лапласа; 

е реальна поверхня Землі не співпадає з референц-еліпсоїдом, тому гранична умова в задачі Молоденського наближена, 1 нев'язка, 
яка містить інформацію про похибку задачі Молоденського для заданої моделі, неліквідна, і сягає - 7 мГал, або понад 290 від пов- 
ної спостереженої аномалії сили тяжіння (ця задача генерує досить грубі наближення); 

е не існує способу вимірів гравітаційного потенціалу, а виміри його похідних на поверхні Землі дорогі, зате накопичені у великих 
обсягах дані гравімагнітних зйомок, що значеннями приростів модуля градієнта потенціалу; 

. у класичних задачах ці дані можна використати лише як наближені граничні умови, що забезпечує визначення потенціалу з гаран- 
тованою точністю лише в локальній області малої міри; 

. відновлення потенціалу в глобальній області цим шляхом приречене на невдачу - розв'язки відповідних граничних задач визначені 
в кожній локальній області з точністю до невизначеної сталої, залежної від геометрії (розмірів і форми) локальної області; 

е | для,склеювання" локальних рішень відсутні які-небудь розумні критерії; 

е підгравітаційним полем розуміють вертикальну похідну потенціалу сили тяжіння, тоді як воно тотожне значенням модуля градієн- 
та потенціалу сили тяжіння; 

. напрям вектора сили тяжіння невідомий, через що задачу перерахунку поля не можна в загальному випадку звести до лінійних за- 
дач математичної фізики. 

Тому задача відновлення потенціалу за значеннями модуля його градієнта, сформульована М.А. Алексідзе, є більш, 
ніж актуальною. Алексідзе вважав, що ,щвивчення питань розв'язності та єдиності цієї граничної задачі ... має принципо- 
ве значення для гравіметрії, оскільки відповість на питання: чи можливе точне однозначне обчислення вектора при- 
скорення сили тяжіння в зовнішньому просторі на основі точного визначення модуля вектора прискорення сили тяжіння 
на замкнутій поверхні, яка цілком охоплює всі тяжіючі маси, і яка мінімальна додаткова інформація потрібна для 
однозначності?" |1, с. 65) (слайд 3). Він довів, що у плоскому випадку задача знаходження гармонічної функції, модуль 
градієнта якої задано на замкнутому контурі, не має єдиного розв'язку. 

Нині у нас є можливість позитивно відповісти на питання, поставлене Алексідзе, та внести суттєві доповнення щодо 
нюансів розв'язності задачі Алексідзе та акцентів її практичного застосування. Але перед цим варто детальніше означити 
деякі результати математичного моделювання аномалій сили тяжіння. 

Властивості модуля градієнта потенціалу (МГПСТ) досліджувалися епізодично, 1 не в останню чергу, через те, що до- 
сі панує необгрунтована точка зору, ніби виміряні аномалії сили тяжіння - це граничні значення або гармонічних функцій 
(у задачах гравірозвідки), або гармонічних функцій 1 їхніх похідних (у задачах геодезичної гравіметрії). У задачах граві- 
розвідки, де розглядають області малої міри (щодо розмірів Землі), такий підхід ефективний і доцільний, а в інших зада- 
чах він вимагає спеціального обгрунтування. 

Розглядаючи аномалії сили тяжіння, закартовані на невеликих ділянках земної поверхні і обумовлені неглибоко роз- 
ташованими джерелами, К. Юнг рекомендував вважати їх у першому наближенні значеннями ,вертикальної" складової 


У, притяжіння збурювального тіла (тобто, замінити граничні дані у вигляді значень МГПСТ значеннями гармонічної 


функції потенціалу У, ). При цьому не зауважено, що похибка такої підміни ростиме разом з розмірами області, де зосе- 


реджені аномалієтвірні об'єкти. Таке тлумачення гравіаномалій в епоху становлення геофізичних методів відкрило широкі 
можливості в інтерпретації даних за розвинутою теорією гармонічних 1 аналітичних функцій. З тих пір граничні задачі Ді- 
ріхле й Неймана для рівняння Лапласа зайняли чільне місце в теорії прямих і обернених задач граві- і магніторозвідки. 

Захоплення теорією, збудованою на базі лінійних граничних задач потенціалу 1 гармонічних функцій настільки си- 
льне, що у дослідники лишили поза увагою зауваження М. А. Алексідзе в задачі редукування значень сили тяжіння: ано- 
малії сили тяжіння, що сприймалися всіма від К. Юнга як гармонічні функції, такими по суті не є. 





Поверхня Землі | 





Рис. 1. Розходження векторів реальної і нормальної сили тяжіння. 

Точні граничні умови в класичних задачах відновлення потенціалу формуються слідами потенціалу або його похід- 
них. Але виміри потенціалу або вектора напруженості сили тяжіння на фізичній поверхні Землі пов'язані з надзвичайними 
технічними і економічними труднощами, нездоланними у осяжному майбутньому. Накопичені гравіметричні дані, але за 
допомогою гравіметрів виміряють тільки прирости МГПСТ, а не значення певної складової градієнта потенціалу, як це ча- 
сто розуміють. При вимірах гравіметри встановлюють в кожній точці земної поверхні за рівнем на еквіпотенціальній 
поверхні, що перетинає поверхню Землі в цій точці (слайд 4), але просторова орієнтація гравіметрів ніяк не фіксується і 
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лишаються невідомими нахили приладів при їхньому зміщенні від пункту до пункту на профілі. Нахили приладів, обумов- 
лені різною кривизною кожної з еквіпотенціальних поверхонь, що проходять через зазначені пункти і непаралельних одна 
одній, можна описати через приріст кута між нормалями перетнутих у точці поверхонь Землі і еквіпотенціальної, або ж 
безпосередніми вимірами значень спрямовуючи косинусів. Але ні ці, ні жодні аналогічні величини, які визначають напрям 
сили тяжіння, не виміряють в точках земного рельєфу через виняткову складність таких спостережень в польових умовах. 
Розходження цих напрямів в реальних умовах коливається від 0" до 40" і може призвести до того, що глибинні 
неоднорідності, які зумовлюють ці розходження, не відображатимуться в класичних аномаліях. 

Реальний результат гравіметричних вимірів - прирости абсолютних значень прискорення сили тяжіння, які відповіда- 
ють приросту модуля вектора сили тяжіння між еквіпотенціальними поверхнями у точках земного рельєфу. Якби на по- 
верхні Землі (рівняння якої є заданим) крім значень МГПСТ і внутрішньої нормалі до поверхні вимірювався і напрям гра- 
дієнта, нелінійну граничну задачу відновлення потенціалу сили тяжіння можна було б звести до розв'язання лінійної зов- 
нішньої задачі Неймана для рівняння Лапласа. 

Самі значення МГПСТ для простого шару, як доведено працями Якимчика, негармонічні, неперервні, на ляпуновській 
поверхні задовольняють диференціальному рівнянню сили тяжіння типу Клейна-Гордона (слайд 5) 
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зі змінним коефіцієнтом а" (Х) , що описує кривизну поля. Якби функція а" (Х) була відома, воно стало б лінійним, і для 


нього з граничними даними значень сили тяжіння легко б вирішувалась задача Діріхле для замкнутої області - зручний 
апарат перетворень МГПСТ, що виникають у прикладних задачах теорій фігури Землі 1 тлумачення гравіаномалій. Однак, 
маємо у розпорядженні лише обмежену інформацію про цю функцію, що відразу ж істотно звужує можливості даного 


способу перетворень 2(Х). 


При перетині густинних границь перші і другі похідні МГПСТ зазнають розривів неперервності, а самі значення 
МГПСТ лишаються неперервними, а їхні особливі точки на границях мають вигляд , кутових", ,повернення", , самопере- 
тину", тощо. 

На підставі даних 2(Х) не можна сконструювати точні граничні дані для жодної з лінійних класичних задач теорії 


потенціалу, у тому числі й для задач Стокса-Молоденського для рівняння Лапласа або Діріхле для рівняння сили тяжіння. 
Для точного рівняння в глобальній області можна сформулювати лише наближені граничні умови, або скласти тільки 
наближене лінійне рівняння для точних граничних даних. 

В обох випадках не можна одержати з високою точністю розв'язок задачі в глобальній області; його пошук - першо- 
чергова проблема теорії МГПСТ. Зауважмо особливість означення аномалії в теорії МГПСТ (слайд б). 


Загалом уточнену модель гравіаномалії в теорії МГПСТ описує вираз до(х) - 2(х)- /(90 -2209) віп? жо де У(х) 


у 
- нормальне значення сили тяжіння, яке відповідає відомому розподілу мас усередині сфероїда, МК) - кут між нормаля- 


ми у точці спостережень (відхилення виска). Якщо нормальний сфероїд 1 розподіл мас у ньому обрані вдало, відхилення 
виска у районах зі , спокійним" рельєфом не перевищуває 60", 1 оцінка зверху останнього доданка у означенні аномалії не 


перевищує величини 2(х)зіп? 209 « 0,05 мГал. Тому аномалії задані як різниця д2(Х) -- 2(х) - У(Х), що за виглядом 
2 


співпадає з класичним означенням аномалії, але істотно відрізняється тим, що нормальні значення сили тяжіння тісно зв'я- 
зані з конкретним - - правильним --- розподілом мас всередині сфероїда. 

Це означення аномалії відрізняється від класичного: в класичній схемі аномалії визначають як відхилення реального 
поля сили тяжіння Землі від деякого його осереднення, яке є певною апроксимацією поля сфероїда з невідомим розподілом 
має усередині нього. Відсутність інформації про розподіл нормальних мас ускладнює інтерпретацію гравіаномалій. 

Похибка продовження аномалій сили тяжіння внаслідок заміни субгармонічної функції гармонічною функцією зале- 
жить від кривизни еквіпотенціальних поверхонь поля, величини аномалій 1 міри локальної області, у яку здійснено продов- 
ження (слайд 7). По суті, це заміна рівняння сили тяжіння типу Клейна-Гордона рівнянням Лапласа. За відомої амплітуди 
аномалій похибки відхилення субгармонічної функції від гармонічної функції у локальних областях коливаються в межах 
4.78 « є «290,6 мГал. Похибка продовжених аномалій як гармонічних функцій в зовнішню область за умови, що амплі- 
туди продовжуваних аномалій не перевершують 150 мгал, для сфероїда, який апроксимує земну поверхню, з параметрами 
а - 6378,2 км, Б - 6356,9 км, матимемо 3,67 « є « 73,41мГал. При цьому вказані величини ще можуть бути занижені. 


Порівнюючи точне аналітичне значення аномалії для сфери з обчисленим з інтегралу Пуассона, маємо відносну похибку г 
3390, тому методи аналітичного продовження гармонічних функцій для локалізації аномалій від окремих тіл, незастосовні 
в глобальній області. Отже, похибки заміни функції, що задовольняє в локальній області рівнянню аномалій сили тяжіння, 
функцією, гармонічною в тій же області, сягають, залежно від виду трансформації, недопустимих величин. 

Оцінки похибок трансформацій в локальних областях, що навіть перетинаються між собою, не дають змоги оцінити 


загальну похибку трансформації в об'єднанні областей. Якщо в одній точці відхилення виска (ХХ) -60", а в іншій точці 


Зх) -0", то при переході від однієї локальної області до іншої істинна аномалія відрізнятиметься від класичної на вели- 


чину 5а(х) - 2У(х)8іп ви я. 285,25 мГал, навіть якщо класична аномалія дорівнює нулю в обох точках. 

Перетворення поля в кожній з областей при однаково рівних нулю класичних аномаліях за правилами трансформації 
гармонічних функцій не дає можливості виявити істинну аномалію, яка сягає величин, порівняних зі значеннями трансфо- 
рмованих аномалій навіть у районах зі спокійним рельєфом. Апроксимації аномалій гармонічними функціями з гарантова- 
ною точністю можна виконати тільки в локальних областях, тому методи гравіметрії, попри їхню розмаїтість 1 загальність, 
дають змогу ефективно вирішувати за заданими аномаліями ту чи іншу прикладну задачу винятково в областях малої міри. 
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Локальну область застосування методів інтерпретації аномалій сили тяжіння, розвинених на базі теорії гармонічних 
функцій потенціалу, можна поширити на вивчення глибинної регіональної структури Землі лише шляхом аналітичного 
продовження модуля його градієнта. Через це задача відновлення потенціалу за значеннями модуля його градієнта - одна 
з найактуальніших проблем гравіметрії. 

Кризова ситуація, що виросла з невідповідності можливостей лінійної апроксимації аномалій, не усвідомлена інтер- 
претаторами і теоретиками, які всупереч дослідженням М.А. Алексідзе, завзято не бачать ,білі плями" у використанні 
методів теорії лінійних граничних задач потенціалу для тлумачення аномалій МГПСТ. Нелінійна гранична задача визна- 
чення потенціалу за значеннями МГПСТ, позитивне розв'язання якої зняло б усі накопичені проблеми, досі не посіла на- 
лежного їй місця серед проблем теорії інтерпретації геофізичних вимірів. 

Уперше її сформулював М. А. Алексідзе як для зовнішню нелінійну задачу для точного рівняння Лапласа 1 наближе- 
них граничних даних - виміряних значень сили тяжіння (МГПСТ) на обмеженій ділянці земного рельєфу, що є частиною 
границі локальної області, для якої ставилася задача, і ,вгаданих" значень сили тяжіння на інших границях, які об- 
межують область (слайд 8). В такому вигляді вона не має єдиного розв'язку. 

Для успішного її розв'язання А.В. Чорний запропонував таку альтернативу (слапіо 9): 

е відшукання такого точного диференціального оператора, який анулює значення модуля градієнта потенціалу поза областю розта- 
шування тяжіючих мас, і розв'язання для цього оператора відповідної лінійної граничної задачі для рівняння Лапласа; 
е постановка такої нелінійної граничної задачі для рівняння Лапласа, в крайових умовах якої безпосередньо задіяні точні значення 

СИЛИ ТЯЖІННЯ. 

Результатом першого підходу є постановка 1 розв'язання зовнішньої задачі Діріхле з даними на поверхні Ляпунова 
для лінійного диференціального рівняння типу Клейна-Гордона (яке моделює значення аномалій МГПСТ в області, не зай- 
нятій тяжіючими масами) з невідомим змінним коефіцієнтом, що відповідає обраному нормальному потенціалу. Цей наб- 
лижений спосіб (бо змінний коефіцієнт рівняння сили тяжіння, для якого збудована теорія, відомий лише першими набли- 
женнями, які характеризують кривизни нормальних полів з еквіпотенціальними поверхнями у вигляді сімейств концент- 
ричних сфер і конфокальних еліпсоїдів) модернізований у дисертації А.І. Якимчика. Для відновлення потенціалу притя- 
жіння за значеннями МГПСТ за умови, що відновлюваний потенціал не надто ухиляється від заданого, знайдена послідов- 
ність розв'язків задачі Неймана для рівняння Лапласа, які визначають збурювальний потенціал. 

Інший підхід якраз і полягає у вирішенні сформульованої Алексідзе задачі, але з граничними даними на класі повер- 
хонь Ляпунова за умови, що відновлюваний потенціал не дуже відхиляється від заданого. Доведено, що задача Алексідзе 
для рівняння Лапласа коректна на множинах витокоподібно поданих функцій потенціалів простого шару, розповсюдже- 
них на випуклих поверхнях Ляпунова (слато 10). 

Для визначення густини потенціалу простого шару, за допомогою якого однозначно розв'язується задача Алексідзе 
розробено ряд алгоритмів (слайд І Г) 


Алгоритм 1. сь (х)з Гкуьбоа, ус» (25. зх с), 
рю 


оскільки ядро рівняння 1 його права частина залежать від густини, для її визначення придатний ітераційний процес 


Со (х) з з(х), С»д (2) -0,хеду, бін (х) 2-6), (х)чн Си (У), ори (х) -А Ід», (9) п 0,со 
де А, РЕЯИ (2) і Бо Си )- Їк, (» енер у», (2)45, 


ду 


Алгоритм 2. ст, (х у Ї к ( х, бус )сь (2)15, є Б (ж, су)» де В, (зо)з (е202- Е(хо)- Ес» Уу с (х). 
у 
для визначення послідовних наближень густини потенціалу (3) можна запропонувати такий ітераційний процес 


10 (2) - з(х) ерже (х) -0,хєду, са (х)- о (хе ери (9), дані (х) -4, Ід», (2) побоо, 
де ДА, дол (9 зб (а Ол )- Гк, (о роди у», (245, 


ду 


Алгоритм 3. Флн(х)З Гк, (х, пан о У 5. зб (а сі, бі (х) золи ()- є, (х), па0,со; 


бу 
оо (х) - 2(х), Од (х) - 0, хє бу . 
Практичні застосування розв'язків задачі Алексідзе можуть бути наступними (слайд 12) 


Уточнення фігури Землі. На поверхні Землі задані значення 7 (х) і початкове наближення фігури Землі О".на і-му наближен- 





ні поверхні о! збурювальний потенціал дорівнює УЧазхьЇ- Т 4 у з. 05 (х), х є о», х, є бу» де (/, (х) - 
ло уг Й 


нормальний потенціал, густина О; (х, , Х» Її; х є 9", х, є ду є розв'язком рівняння 


а усбанго) сов(т, р), (2,х,)є епьх ни і 
аа з я Ї : п хз 8, нЖЦуз Ії ри я (р, аа8 48, з 2"(,) 


а аномальна висота за відомою формулою Брунса о (х ЧЕ аое У (о // УС), х Є о", Х, є ду визначає поверхню 11 /-го на- 








ду оз Ї" дуду 


ближення форми Землі за формулою о0"'зб'з ло" Є і . На і-му кроці обчислюємо О" - наближення поверхні Землі, густину 
С; о і потенціал У"(хух), за формулою Брунса - наближення форми і ре З Х є б" геоїда, згодом - наступне 


6" Е - і и 
наближення 1 так далі - до заданої точності. 


Відновлення магнітного потенціалу. Якщо магнітний потенціал - сума УС)-п (х)-- Т(х) нормального С (х) і 
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збурювального т(х) магнітного потенціалу для даної епохи, відновлення магнітного потенціалу У(2) о За заданими на границі 
Хан 


, можна звести до ітераційного процесу відновлення збурювального потенціалу 





значеннями модуля його градієнта 7 (х) ш | ггад У(Х) 


9 2 
внусу ОО), дТ,(), ук цу - 902) , д"Т, (0) 
М А хо дж х"(9 дх,дх, | дх,дх, 


й в ба 8 ре А б (ру 


дх, 4х 3 |х- єї дхудх, (4я2 00 ре-еЇ 


9; 

а спрямовуючі косинуси для наступного наближення модуля градієнта магнітного потенціалу 7. , ( х) 27 () -а ) (п; паб ) а 
їч іч іч 

ко 


, 














іні - ; : 
де сов(п, пре ) б уб 1/ УДІЕ (х) , потрібен для обчислення чергового наближення збурювального магнітного потенціалу 
Т. а (2) 2, (х)сок(ли т) Ух) собу, т). 
: : (іа) б 

Наступне наближення магнітного потенціалу И (х) -(/ (х) з Т, з (х) и 

Нормальний потенціал - аномальне поле однорідно намагніченої сфери. Окрім даних 7 (х) (АТ, ( х)) відомі рівняння поверхні 
Землі і вектор одиничної нормалі т(х) -- сов(х, т, ) К - 1,2,3 до неї майже в будь-якій її точці, і нормальне поле (/ о 

1д/ (х) 


Нормаль у(х) ш со5(у, Хі ) до еквіпотенціальної поверхні (/ (х) з Сх дорівнює со5(у, Хі )- --3-Х-2 тазначення 
у бх, 


ЕІ 
1 
со5(у, т) - У сову, х, |соз(х,, т), де - 2 Івгад М (2. - 
КІ у (») 
: Зоо . - з ; 
Якщо поверхня Землі є множиною Ляпунова, то магнітний потенціал У(х), Хх є у для даної епохи відновлюється однозначно 
за даними вимірювання величини напруженості 7 (х), х є бу магнітного поля цієї епохи. 


Трансформація сили тяжіння. Для розв'язання задачі трансформації і перерахунку сили тяжіння формулюється така нелінійна 
гранична задача для потенціалу сили тяжіння: 








Для розв'язання цієї задачі методом розкладу за неортогональними функціями застосовано ітераційний процес 


А -0, аукт)07ОЬ),, а (уз) 007Сь У) 5 
дх бу 


х , . ние 
де (/" - К -те наближення розв'язку задачі, а коефіцієнти обчислюються з формул 














2! 2 
Р кад. хх - , ()- і 
о! опе! дО! оо! 
рів АР ОР рев тн, в 
дх бу дх бу 


Результати чисельних експериментів з ітераційного розв'язання нелінійної граничної задачі Алексідзе у плоскому 
випадку (слайд 1 1) були спрямовані саме на перевірку збіжності цього ітераційного процесу для моделі кулі. 
Вирішено такі тестові зовнішні граничні задачі для рівняння Лапласа: 
г 
ЛИ -0, нтаа ОЇ, - г АЦ «0, 


хзу'за 


г г 
аці - - і 
й ; (хо-тУ вубча (хат) вубжа 














; З і ї : Л 2 2 2 
де 5 - горизонтальна площина, г - 10 і 2 - 5 км. Точний розв'язок першої граничної задачі - г/ хо У ОС) ,другазадача точ- 
ного розв'язку не має. 


; -6 б є . . . о одбо з 5 
Розв'язок з точністю 10 першої граничної задачі в 30 точках поверхні 9 на рівномірній мережі з кроком 10 км в області: -25 
км с у 2 25 км, -25 км « х с 15 км досягнуто за 5 ітерацій. Перше наближення для спрямовуючих косинусів задано таким чином: 


Г2Я (0, )- аз (0,)-0, 3 (о, )- 1, і ос (2, )-0.8, а (0, )-0, аз (0 )- 0.6. На практиці точність визначення похідних (/ , 


у. ї С. буде нижчою. Самі градієнти У співпадають з меншою точністю через погану обумовленість задач з граничними умовами, 


з 
що містять похідні (для них несправедливий принцип максимуму). Цю обставину, непов'язану зі збіжністю ітераційного процесу, 
треба враховувати при вирішенні практичних задач. 

Виявлено, що зовнішня гранична задача вирішується із задовільною точністю, якщо напрям градієнта на всій границі 
області складає з внутрішньою нормаллю гострий кут; крім того, похибка розв'язку зумовлена не поганою збіжністю 
ітерацій (їх число не перевищує 10), а поганою обумовленістю взагалі граничної задачі з косою похідною в граничних 
умовах задачі (напрям зовнішньої нормалі не має бути ортогональним до напряму косої похідної, що не виконується для 
всієї нескінченної поверхні у випадку точкового джерела). 

При розв'язанні практичних задач шукаємо не потенціал У/ сили тяжіння, а потенціал Ї/ аномалії сили тяжіння, зв'я- 


заний з И/ рівністю Й - / - У, де У - відомий потенціал нормального значення сили тяжіння. За лінійних наближень 
граничних задач гравіметрії і магнітометрії такий перехід не змінює не лише рівняння, а й вигляд оператора граничних 
умов, який ще потребує додаткових досліджень 

Основні висновки (слайд 13) на даному етапі можна сформулювати так: 

можна вважати закритим питання можливості отримання поля у зовнішньому просторі з як завгодно високою точністю за відомим 
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розподілом значень модуля градієнта потенціалу сили тяжіння. Ця основна нелінійна задача граві- і магнітометрії обгрунтована в 
алгоритмічних схем її розв'язку. Це важливо тому, що нелінійна апроксимація гравіполя не відображена в курсах гравіметрії, а ми 
демонструємо, що вона справедлива для вевного класу областей. Оскільки це зовнішня гранична задача, для її вирішення недоцільно 


вживати універсальні методи розв'язання граничних задач - скінченних різниць чи скінченних елементів. Труднощі чисельного розв'я- 
зання лавиноподібно зростають при складному рельєфі земної поверхні... 
Дякую за увагу! 


